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Quadrature by -Gauss' Method (DOI or SADQl) 

Closed 

18 + 2n 

1 word at (may also be used by auxiliary 

subroutine) 

n(l.8 + T) milliseconds, where T is the 

duration of the auxiliary subroutine. 

See Note. 

When this routine is located at q entry is 

made by 

— sF 

50 pF 



p + 1 



26 qF 



DESCRIPTION 



where s is the location of the auxiliary 
subroutine which computes values of f(u.). 
This routine uses the method of Gauss 
for approximating the integrals 

-1 



n 



Jf(u)du* I a. f(u.) 
i=l x x 



(1) 



The values f(u.) are obtained from an 
auxiliary subroutine which must be prepared 
by the user. The auxiliary subroutine is 

entered in the standard manner with u. in A, 

l * 

and it should place f(u.) in A. The values 

of the u. are in the range < u. < 1. Loca- 
1 l 

tion is used by this routine to store the 
new value of f (u. ) each time control is 
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returned to this routine by the auxiliary- 
routine. 

At the end of the library copy of this 
routine are sections of tape containing 

constants and the values of u. and a. for 

1 1 

each value of n from 2 through l6. The 
choice of n is made by copying the appro- 
priate section of tape after copying the 
part containing the basic routine. 
ACCURACY Quadrature is carried out numerically by 

approximating the integral with a finite 
summation 



g(x)dx S L ■ E a j[ f(u ± ) (2) 



^b n 

Ja i=l 

where 



L = b - a 
u = (x - a)/L 
f (u) = g(x) = g(Lu + a) 

Gauss* Method of quadrature consists 

essentially of making an optimum choice of 

the values of both the u. and the a. . When 

l l 

the accuracy "with which values of g(x) can 
be obtained is independent of x, the method 
of Gauss gives the highest possible order of 
accuracy for a given number n of ordinates. 
When n ordinates are used, the order of 
accuracy achieved is 2n - 1 $ i.e* the 
truncation error is zero for polynomials of 
degree 2n - 1 or less . For an arbitrary 
function having a continous 2n derivative the 



• 3 - 



truncation error is, for the integral 
shown in equation (l), 



Kfr) 



.(2n) 



ill 
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(nOl 
(2n)l 



n * ' (2n + 1)1 

where 5 is some point in [0,1], and for the 
integral shown in equation (2) 



. / x ,. ,2n+l g (2n) (^) 
R (g) = (b - a) 2 Si — L 

(2n + 1)1 
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l2n7T 



where "^ is some point in [a,b]. 

In addition to the truncation error 
there is a round-off error which may be 
significant. It is given by (see equation (l)) 



= \ [ a.Af(u.) ± 2" lK) f(u.) ± 2^ a.f u.)] + 2-" 4 ° 
i=l X 



l — 
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EXAMPLE 



< +2~ k0 tfax (n|f(u.)|+|f r (u.)|) + A + 2'^ 
(Ku<L X 1 

where A = Max Af (u) is the accuracy of the 

~k0 
auxiliary subroutine. The terms 2 f (u. ) 

~k-Q 1 

and +2 a.f* (u. ) represent errors resulting 

from a maximum error of 2 in the value of 

each of the a. and u., respectively. The 

-kQ 
last term, 2 , represents the final round- 
off of the double precision sum of products. 
The total error is bounded by the sum of the 
truncation error and the round -off error. 
Using n = 5 ordinates in approximating the 



- k - 



integral 
■1 



l 



1 • 1 

s- sin itxdx = — 

2 * 



this routine gave a value of 

0.31850990373, which is off 1 unit 

in the 8 decimal place. 

Using the same number of ordinates integration 

"by Simpson's rule (code E-3) gave a value of 

0.31903559373, which is off 1 unit in the 

third decimal place. 
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00 K(E5) 











L5 i6l 
46 % 




Set addresses 
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42 JL 
K5 F 
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46 Gl 

42 13L 




Plant S 
Plant link 
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41 i4l 
2S 4l 




DSt £-1 = C. 
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S5 F 
40 15L 


from 12 
from 3 
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L5 F 
50 5L 


by 0_, 10 


'"1 

Enter aux* routine 
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26 F 
40 F 


by 2 


f(u ± ) 
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L5 15L 
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50 F 
74 F 
L4 i4l 


by },>■ 11 


E a.f (u.) 

. , 11 
i=l 




9 


4o i4l 
L5 8L 




Step addresses 
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L4 5L 
46 5L 
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F5 TL 
^2 TL 
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LO 17L 
36 4L 




Test i > N 
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J f(u)du 
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by 2 


link 
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00 F 
00 F 
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00 F 
00 F 
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E 5 
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1-7 
18 

19 



n + 17 
n + 18 
n + 19 



00 18L 

00 (n+l8)L 

75 15L 

50 (2n + 18 )L 

U l 



n 



2n + 17 



